Quantum walks of correlated particles 
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Quantum walks of correlated particles offer the possibility to study large-scale quantum interfer- 
ence, simulate biological, chemical and physical systems, and a route to universal quantum compu- 
tation. Here we demonstrate quantum walks of two identical photons in an array of 21 continuously 
evanescent ly-coupled waveguides in a SiOxN^ chip. We observe quantum correlations, violating a 
classical limit by 76 standard deviations, and find that they depend critically on the input state 
of the quantum walk. These results open the way to a powerful approach to quantum walks using 
correlated particles to encode information in an exponentially larger state space. 



With origins dating back to observations by Lucretius 
in 60BC and Brown in the 1800's, random walks are a 
powerful tool used in a broad range of fields from genetics 
to economics |jLj. The quantum mechanical analogue — 
quantum walks [2j |3] — corresponds to the tunnelling of 
quantum particles into several possible sites, generating 
large coherent superposition states and allowing mas- 
sive parallelism in exploring multiple trajectories through 
a given connected graph {eg. Fig. [T]). This quantum 
state evolution is a reversible (unitary) process and so 
requires low noise (decoherence) systems for observa- 
tion. In contrast to the diffusive behaviour of (classical) 
random walks, which tend towards a steady state, the 
wave function in a quantum walk propagates ballistically 
(Fig. |2|c)). These features are at the heart of new algo- 
rithms for database-search [4^ , random graph navigation, 
models for quantum communication using spin chains [5] , 
universal quantum computation [6 and quantum simu- 
lation [7 . 

Quantum walks have been demonstrated using nuclear 
magnetic resonance [H [9] , phase [TOl |TT] and position [12] 
space of trapped ions, the frequency space of an optical 
resonator [13 , single photons in bulk [14 and fibre 
optics and the scattering of light in coupled waveguide ar- 
rays ^6|. However, to date, all realisations have been lim- 
ited to single particle quantum walks, which have an ex- 
act mapping to classical wave phenomena [17 , and there- 
fore cannot provide any advantage from quantum effects 
(note that the quantum walk with two trapped ions [11] 
encodes in the centre of mass mode and is therefore effec- 
tively a single particle quantum walk on a line). Indeed 
single particle quantum walks have been observed using 
classical light [TgIIT^. In contrast, for quantum walks of 
more than one indistinguishable particle, classical theory 
no longer provides a sufficient description — quantum the- 
ory predicts that probability amplitudes interfere lead- 
ing to distinctly non-classical correlations [19^ i20| . This 



quantum behaviour gives rise to a computational advan- 
tage in quantum walks of two identical particles, which 
can be used to solve the graph isomorphism problem for 
example [21]. The major challenge associated with realis- 
ing quantum walks of correlated particles is the need for a 
low decoherence system that preserves their non-classical 
features. 

The intrinsically low decoherence properties and easy 
manipulation of single photons make them ideal for ob- 
serving quantum mechanical behaviour and for quan- 
tum technologies [22]; and the effectiveness of arrays 
of continuously coupled waveguides for bright classical 
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FIG. 1: Quantum walks with one and two indistinguishable pho- 
tons described by a Hamiltonian of coupled harmonic oscillators 
(Eq. [ij). (a). The linear array of vertices representing the state 
space of one photon populating N = 21 waveguides in a coupled 
one-dimensional array. Site potentials for each vertex are — /3 while 
hopping amplitudes between nearest neighbours are all equal to 
—C. (b). The two-dimensional lattice of vertices that represent 
the state space of two photons populating N = 21 waveguides in a 
coupled one-dimensional waveguide array, (inset) Enlarged portion 
of the lattice displays the vertex representation of the two photon 
basis states. Site potentials are all equal to —2/5, while hopping 
amplitudes between adjacent vertices are either —C or —^/2C as 
labeled. 
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light has been demonstrated [161 [18]. However, quantum 
walks with correlated photons in such structures requires 
a means to measure two-photon correlations across the 
waveguide array. The spacing between waveguides in an 
array required for evanescent coupling (of order several 
/im) is significantly smaller than the minimum spacing 
of optical fibre arrays (127 fim) typically used to cou- 
ple to single photon detectors. Previous quantum opti- 
cal waveguide circuits in a silica-on-silicon architecture 
[22 promise to avoid decoherence effects [23] present in 
other experimental realisations, interferometric stability 
and near-perfect mode overlap [22], which is problem- 
atic in large scale bulk optical realisations (an arbitrary 
N mode multi-port would require a network of 0{N'^) 
beams splitters [24 ). However, the low refractive in- 
dex contrast (A = (n^^^^ - n%^^^.^g)/2nl^^^ ^ 0.5%) in 
this architecture results in a large minimum bend radius 
(< O.ldB loss at 800 nm) of ^ 15 mm, making it un- 
suitable for coupling into and out of large array quantum 
walk devices. 

We have overcome these technical challenges by fabri- 
cating waveguide arrays in SiO^gN^ (silicon oxynitride), 
a material that enables a much higher refractive index 
contrast than silica-on-silicon (the refractive index is de- 
termined by X and y), resulting in more compact de- 
vices (Fig. [2|a)) and a practical means to realise large 
coupled waveguide arrays that can be coupled to optical 
fibres. The device shown in Fig. [2|a) is a 5 mm long 
silicon chip with SiOa^N^y waveguides with high refractive 
index contrast A = 4.4%. The minimum bend radius for 
this index contrast is 600 /am which enables much more 
rapid spreading of the waveguides from the evanescent 
coupling region where waveguides are pitched at 2.8/im 
to a pitch suitable for optical fibre (250/im and 125/im 
for photon injection and photon collection respectively, 
see Appendix). 

Photon pairs were generated via type-I SPDC in a 2 
mm thick nonlinear bismuth borate BiBaOe (BiBO) 
crystal, pumped with 40 mW of 402 nm light from a 
CW diode laser. With an opening angle of 3°, pairs of 
degenerate A = 804 nm photons were filtered by 2 nm 
FWHM interference filters, and focused onto two inputs 
of an array of polarization maintaining fibres which were 
butt-coupled to the waveguide chip. At the output of 
the chip, an array of multi mode fibres guide the out- 
put to 12 single photon counting modules connected to 
photon counting logic based on three FPGA boards used 
to measure matrices of two-photon correlations between 
the outputs of the array (Figs, [s] and |4|. Quantum inter- 
ference (degree of indistinguishability) was controlled by 
relative temporal delay between the pair of photons, us- 
ing an automated linear actuator, and characterised with 
a standard Hong-Ou-Mandel experiment ^25 . The pitch 
of the waveguide outputs is half of that of the collecting 
fibre arrays allowing 121 of the 231 possible two-photon 
correlations to be measured at the output. 
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FIG. 2: A continuously coupled waveguide array for realising 
correlated photon quantum walks, (a) An optical micrograph 
of a 21 waveguide array showing the three input waveguides, 
initially separated by 250/im bending into the 700/im long 
coupling region. All 21 outputs bend out to 125/im spacing, 
(b) Simulation of single photon propagation in the array, (c) 
Output pattern of 810nm laser light propagating through the 
waveguide array. 

Photons propagating through the coupled waveguide 
array shown in Fig. 2(a) are modelled assuming nearest 
neighbour interaction with the Hamiltonian for coupled 
oscillators [26 , setting h = 1: 



H 



N 



where the creation and annihilation operators Oj and aj 
obey Bose-Einstein statistics and act on waveguide j. In 
our devices the waveguide propagation constant Pj = 
/3 and coupling constant between adjacent waveguides 
^jj±i = C = C are designed to be uniform for all j. 
Through choice of operator A the Heisenberg equation of 
motion idA/dz = [A, H] models the dynamics of photons 
propagating along distance z of the array. For example, 
injecting single photons in waveguide k is described using 



(1) 



yielding [27j 
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Acting Eq. [2] on the vacuum |0), and inspecting in the 
Schrodinger picture idltpi) /dz = H'^^^ yields H'^^^ 
From this Hamiltonian, an adjacency matrix Mj'^jj. = 
j (1| H^^'^ |1)^ of the graph in Fig. ma) is constructed us- 
ing the single photon Fock basis {|l/j} representation. 
A quantum walk on this graph evolves over length z ac- 
cording to the unitary transform /7^^^ = exp[— zi^^^^z] 
[31 [27], the exact dynamics of which can be observed 
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FIG. 3: Measured and simulated correlations in waveguide arrays when two photons are coupled to waveguides and 1: 
ajaj |0) (a,b) for input photons separated with temporal delay longer than their coherence length and (c,d) for photons arriving 
simultaneously in the array. All resulting measurements are corrected for coupling fluctuations using simultaneously detected 
single photon signal and as well relative detector efficiency; the integration time was 1 hour. The outcome of two photons 
populating one waveguide was detected using non-deterministic photon number resolving detection using an optical fibre splitter. 



by injecting bright light into the waveguide array [16] . 
We used this approach to calibrate our device and the 
coupling efficiencies by launching horizontally polarised 
810 nm laser light into the input of the central waveg- 
uide, waveguide 0. We measured the interference pattern 
shown in Fig. [2|c) which corresponds to the probability 
distribution for single photons detected at the 21 output 
waveguides (numbered -10 on the left, through 0, to 10 on 
the right). Before the waveguides reach their minimum 
separation at the central coupling region, they are sig- 
nificantly coupled in the spreading regions on both sides. 
This coupling provides an effective coupling length of 82 
jam. We determined this length by comparing the output 
interference pattern (not shown) for an array of length 
350 /im length. From these data the coupling constant 
C = 5 mm~^ was determined by running the simulation 
and minimising the square of the errors. 

The injection of two indistinguishable photons is mod- 



elled with the operator A 



dajat , , 

dz 3 



a]al_ 



a^-a\ yielding: 



1+44+1+44-1+44+1 



(3) 



From Eq. [sj H^'^^ acting on the two photon Fock space 
is extracted as a matrix represented in the two-photon 
Fock basis {|1)^- |1)^ , |2)^} which is equal to the adja- 
cency matrix for the graph given in Fig. [ijb). Evolu- 
tion of the two-photon state through the device therefore 
simulates a single particle quantum walk on the graph 
in Fig. [ijb) with 0{N'^) vertices and unitary transform 
= exp[— zi^^^^z]. In general, linear increase in the 
number of photons input into the coupled array results 
in exponential growth of the Hilbert space and of the 
corresponding graph on which a single particle quantum 
walk is emulated. However, only when indistinguishable 
photons are injected in the device can the output state 
be non-seperable, exhibiting non-classical correlated be- 
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FIG. 4: Measured and simulated correlations in waveguide arrays when two photons are coupled to waveguides -1 and 1: 
aLittj |0) (a,b) for input photons separated with temporal delay longer than their coherence length and (c,d) for the photons 
arriving simultaneously in the array. All measurements are corrected for coupling fluctuations using the simultaneously detected 
single photon count rate and relative detector efficiency; the integration time was 1 hour. The outcome of two photons populating 
one waveguide was detected using non-deterministic photon number resolving detection using an optical fibre splitter. 



haviour (see Appendix). Note that the two photon uni- 
tary evolution can also be computed from the product of 
single photon mode transformations [27] . 

The measured correlation matrices F^^^ (defined as 
probability to detect a two photon coincidence across 
waveguides r) ^1 ^ for injecting two single pho- 
tons into the central neighbouring waveguides and 1 
{i.e. ajaj |0)) are plotted in Fig. [sj^a) for photons 
made distinguishable using temporal delay (not over- 
lapped) and in Fig. [sjc) for pairs of indistinguish- 
able (overlapped) photons. The overlap of these mea- 
sured distributions with ideal simulations (plotted in 
Fig. Isfb) and ^d)) are S = 0.980 ± 0.001 and S = 
0.934 ±0.001, respectively, where S is the similarity be- 
tween two probability distributions F and F' defined by 

S = (Ei,, v''f~r^)VEi,, Ti,, E.J which is a gen- 
eralisation of the average fidelity based on the (classical) 



fidelity between probability distributions. The lower S in 
the overlapped case is attributed to imperfect quantum 
interference. These results clearly display a generalised 
bunching behaviour (tending to both travel to one side of 
the array or the other), characteristic of quantum inter- 
ference: the vanishing of the two off-diagonal lobes is a 
result of destructive interference of quantum amplitudes 
resulting from repeated 7r/2 phase shifts in the photon 
tunnelling between neighbouring waveguides. 

Distinctly different behaviour is observed on injecting 
two photons in two waveguides with one waveguide sep- 
arating them. The measured correlation matrices for in- 
jecting photons into waveguides —1 and 1, and the vac- 
uum in waveguide {i.e. aLi«l |0)) in the centre of the 
array are plotted in Fig[4j The similarities with the ideal 
simulation are S = 0.970 ± 0.002 and S = 0.903 ± 0.002 
for the delayed and overlapped photons, respectively. In 
this case, instead of bunching, when the two photons are 



FIG. 5: Violating the classical limit. Number of standard deviations a of the violation of inequality Eq. [5] for injecting 
indistinguishable photons into the inputs 0, 1 (a) and -1, 1 (b) and measuring correlations at waveguide outputs labelled q 
and r. 



indistinguishable they generate a pattern where the main 
feature is the vanishing of probability to simultaneously 
detect one photon in the centre of the array and a one at 
the limit of ballistic propagation (for example in waveg- 
uides and 7 in Fig|4jc),(d)). 

Detecting two-fold coincidences of two indistinguish- 
able photons leads to non-classical correlations across 
pairs of waveguide in the array [27] . The correlation func- 
tion after length z for two photons populating waveguides 
q and r is given [271 HH] by: 



1 



^(1) ^(1) , ^(1) ^(1) 
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For classical light, including random phase fluctuations 
that mimic certain properties of quantum light, diagonal 
correlations Tq^q are related to correlations in the off- 
diagonal lobes Tq^ri q 7^ r according to the inequality 
EH: 



.(CO p(co 



> 0. 



(5) 



Inequality Eq. [5] is violated when two indistinguish- 
able photons are injected into the device. Measured vi- 
olations from injecting photons into 0, 1 and -1,+1 are 
plotted Fig. |5j with white data points representing no vi- 
olation, and coloured data points representing the extent 
of violating Eq. [5] for each pair of waveguides. This is 
quantified as a function of standard deviations a (com- 
puted from propagation of error from two photon coin- 
cidence detection, assuming Poissonian statistics), with 
the maximum violation reaching 76 standard deviations. 
Inequality Eq. |5] is not violated when the photons of the 
input pairs are distinguishable. 

These demonstrations show uniquely non-classical be- 
haviour of two identical particles, tunnelling through ar- 
rayed potential wells; two photons initially prepared in 
a separable product state interfere in a generalisation 
of the Hong-Ou-Mandel effect [25' yielding non-classical 



spatial correlations. Increasing the photon number n will 
emulate quantum walks on hyper-cubic graphs exponen- 
tially large in n, while exploiting quantum interference 
in three-dimensional directly written waveguides [29^, [30] 
allows further increase in graph size [31 . Enlarging the 
guided mode or decreasing the waveguide separation pro- 
vides another increase in graph complexity. This requires 
a model beyond nearest neighbour coupling, tending to- 
wards a multimode interference slab waveguide when the 
channel separation goes to zero [32 . Here, we have mod- 
elled and used arrays of waveguide with fixed propaga- 
tion and coupling coefficients = P and Ck±i = C, 
but varying these parameters independently for each k 
provides a means to engineer the quantum walk's precise 
graph structure. For example, varying these parameters 
randomly and independently allows investigation of cor- 
related quantum walks in disordered systems and verify 
the effects of Anderson localisation |33], known to affect 
propagation of quantum information [34l 13^ . Reconfig- 
urable waveguide circuits [22] will allow real-time control 
of the input state and the graph structure itself, enabling 
for example phase control over entangled "NOON" in- 
put states to yield methods for simulating symmetric and 
anti-symmetric particles undergoing quantum walk [27] . 

We thank N. Brunner, A. Laing and T. Rudolph for help- 
ful discussion. This work was supported by EPSRC, ERC, 
QUANTIP, QIP IRC, lARPA, the Leverhulme Trust, and 
NSQI. J.L.OB, acknowledges a Royal Society Wolfson Merit 
Award. 



Electronic address: "Jeremy.OBrien@Bristol.ac.uk' 
[1] R. Motwani and P. Raghavan, Randomized Algorithms 

(Cambridge University Press, 1995). 
[2] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. 

A 48, 1687 (1993). 
[3] E. Farhi and S. Gutman, Phys. Rev. A 58, 915 (1998). 



6 



[4] A. M. Childs and J. Goldstone, Phys. Rev. A 70, 022314 
(2004). 

[5] S. Bose, Phys. Rev. Lett. 91, 207901 (2003). 

[6] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2008). 

[7] M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru- 

Guzik, J. Chem. Phys. 129, 174106 (2008). 
[8] C. A. Ryan, M. Laforest, J. C. Boilequ, and R. Laflamme, 

Phys. Rev. A 72, 062317 (2005). 
[9] J. Du, H. Li, X. Xu, M. Shi, J. Wu, X. Zhou, and R. Han, 

Phys. Rev. A 67, 042316 (2003). 
[10] H. Schmitz, R. Matjeschk, C. Schneider, J. Glueckert, 

M. Enderlein, T. Huber, and T. Schaetz, Phys. Rev. Lett. 

103, 090504 (2009). 
[11] F. Zahringer, G. Kirchmair, R. Gerritsma, E. Solano, 

R. Blatt, and C. F. Roos, Phys. Rev. Lett. 104, 100503 

(2010). 

[12] M. Karski, L. Forster, J.-M. Choi, A. Steffen, W. Alt, 
D. Meschede, and A. Widera, Science 325, 174 (2009). 

[13] D. Bouwmeester, L Marzoli, W. Karman, G. P. nd Schle- 
ich, and J. P. Woerdman, Phys. Rev. A 61, 013410 
(1999). 

[14] B. Do, M. Stohler, S. Balasubramanian, and D. Elliott, 

JOSA B 22, 499 (2005). 
[15] A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, 

P. J. Mosley, E. Andersson, I. Jex, and C. Silberhorn, 

Phys. Rev. Lett. 104, 050502 (2010). 
[16] H. B. Perets, Y. Lahini, F. Pozzi, M. Sorel, R. Moran- 

dotti, and Y. Silberberg, Phys. Rev. Lett. 100, 170506 

(2008). 

[17] P. L. Knight, E. Roldan, and J. E. Sipe, Phys. Rev. A 

68, 020301 (2003). 
[18] D. N. Christodoulides, F. Lederer, and Y. Silberberg, 

Nature 424, 817 (2003). 
[19] Y. Omar, N. Paunkovic, L. Sheridan, and S. Bose, Phys. 

Rev. A 74 (2006). 
[20] P. K. Pathak and G. S. Agarwal, Phys. Rev. A 75, 032351 

(2007). 

[21] J. K. Gamble, M. Friesen, D. Zhou, R. Joynt, and S. N. 
Coppersmith, Phys. Rev. A 81 (2010). 

[22] J. L. O'Brien, A. Furusawa, and J. Vuckovic, Nat. Pho- 
tonics 3, 687 (2009). 

[23] V. Kendon, Mathematical Structures in Comp. Sci. 17, 
1169 (2007). 

[24] M. Reck, A. Zeilinger, H. J. Bernstein, and P. Bertani, 
Phys. Rev. Lett. 73, 58 (1994). 

[25] C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 
59, 2044 (1987). 

[26] L. E. Estes, T. H. Keil, and L. M. Narducci, Phys. Rev. 
175, 286 (1968). 

[27] Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silber- 
berg, Phys. Rev. Lett. 102, 253904 (2009). 

[28] K. Mattle, M. Michler, H. Weinfurter, A. Zeilinger, and 
M. Zukowski, Appl. Phys. B 60, Sill (1995). 

[29] G. D. Marshah, A. Politi, J. C. F. Matthews, P. Dekker, 
M. Ams, M. J. Withford, and J. L. O'Brien, Opt. Express 
17, 12546 (2009). 

[30] R. Keil, A. Szameit, F. Dreisow, M. Heinrich, S. Nolte, 
and Tiinnermann, Phys. Rev. A 81, 023834 (2010). 

[31] T. D. Mackay, S. D. Bartlett, L. T. Stephenson, and B. C. 
Sanders, Journal of Physics A: Mathematical and Gen- 
eral 35, 2745 (2002). 

[32] A. Peruzzo, A. Laing, A. Politi, T. Rudolph, and J. L. 
O'Brien, In preperation (2010). 

[33] Y. Lahini, Y. Bromberg, D. N. Christodoulides, and 



Y. Silberberg, arXiv:1003.3657vl (2010). 

[34] J. P. Keating, N. Linden, J. C. F. Matthews, and A. Win- 
ter, Phys. Rev. A 76 (2007). 

[35] Y. Lahini, A. Avidan, F. Pozzi, M. Sorel, R. Moran- 
dotti, D. N. Christodoulides, and Y. Silberberg, Phys. 
Rev. Lett. 100 (2008). 



7 



APPENDIX 

Devices: Each layer defining ttie piiotonic circuit 
was grown on thermally oxidised (8 /im) Silicon (Si 
(100)) wafers with deposition performed in an Oxford 
parallel plate PECVD reactor utilising SiH^ and N2O 
precursors. After annealing channel waveguides are 
obtained by standard lithography and reactive ion 
etching in CHFs/02 chemistry. The channel structures 
are covered by a PECVD silicon oxide {Si02) cladding 
layer. The core and cladding have a refractive index 
contrast of A = (n^^^^ - nli^^^.^g)/2nl^^^ = 4.4% 
enabling a waveguide bend radius of 600 /im resulting 
in compact 5mm long chips shown in Fig. |2|a) of the 
main text. The distance between waveguides in the 
coupled array (2.8 /im) defines the tunnelling rate of 
light between adjacent waveguides, while the waveguides 
at the end of the array spread out, at an equal rate 
of relative separation between adjacent waveguides, to 
a pitch 125 /im for coupling to optical fibre external 
to the chip; each pair of adjacent waveguides bend at 
the same rate until evanescent coupling is negligible, 
ensuring uniform coupling throughout the structure. 
The overall coupling efficiency through the chip is ^10% 
attributed to mode-mismatch between waveguides and 
input/output optical fibre. Injecting single photon input 
from SPDC is achieved using arrays of polarisation 
maintaining fibre, while collection uses arrays of multi- 
mode fibre. Device characterisation used 810 nm light 
from a diode laser source launched into the central input 
with power intensity measured with a CCD camera; the 
resulting interference pattern is given in Fig. ^c) and 
displays the ballistic dynamics equivalent to a single 
photon undergoing a continuous time quantum walk [16] . 

Violations: Violation of the classical limit Eq. |5] 
of the main text is quantified in Fig. |5] by plotting the 
number of standard deviations measured violations reach 
beyond the classical limit. This violation is defined to be 
when the left hand side of Eq. [5]is found to be negative: 

2 / 

^q,r — '^\/^q,qi^r,r < 0- (6) 



The measured value of this quantity for inputs %a^_^-^ |0) 

and aLi^^^L |0) given in figure Fig. ^ (a),(c), from 
which the number of standard deviations of Fig. [5] 
are computed by propagating the Possonian standard 
deviation of measured two-photon statistics. The 
corresponding theoretical model of these violations are 
given in Fig. |6] (b),(d). 

Dimension of mult i- particle quantum walk 
Hilbert space: In general Hilbert space representing n 
distinguishable photons injected into a waveguide array 
of size N is the tensor product of n Hilbert spaces each 
of dimension N. This makes the combined state space 
dimensional. However, the photons still remain in 
a separable product state throughout evolution with 
unitary transform in the form of a tensor product of 
n single photon unitary operators. These dynamics 
are trivial to reconstruct, requiring n uses of the single 
photon evolution (equivalently bright laser light) for a 
uniform array, and n separate evolutions in general. 

For n indistinguishable photons, the dimension of the 
Hilbert space also grows exponentially — bounded below 
by N^/n. However, unlike the distinguishable case, evo- 
lution moves into non-separable states via non-classical 
interference — for example it is well-known that two- 
photon NOON states are generated by non-classical in- 
terference after a 50:50 beam splitter. For the case 
of two indistinguishable photons {A = a^a^) injected 
to array of size the Hilbert space is of dimension 
{N + l)A^/2 (proven by induction). Similarly for 3- 
photons {A = Ojolaj), the Hilbert space is of dimen- 
sion (TV + 2) (TV + l)N/6\ (proven by induction); this 
emulates a CTQW on a three dimensional lattice. For 
n indistinguishable photons non-classically interfering in 
a quantum walk, the dynamics cannot be decomposed 
into the dynamics of n independent walks. The combi- 
nation of distinctly non-classical states created through 
unitary evolution, and the dimension of the Hilbert space 
being exponentially large in n, implies simulation with 
non-interfering photons requires an exponential number 
of resources. 
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FIG. 6: Amount of violations of inequality Eq. [5|of the main text. (a,b): measured and predicted for input state %a^_^-^ |0). 
(c,d): measured and predicted for input state aLia+i |0). 



